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We consider the problem of growing multiplex networks with intrinsic fitness and 
inter-layer coupling. The model comprises two layers; one that incorporates fitness 
and another in which attachments are preferential. In the first layer, attachment 
probabilities are proportional to fitness values, and in the second layer, proportional 
to the sum of degrees in both layers. We provide analytical closed-form solutions for 
the joint distributions of fitness and degrees. We also derive closed-form expressions 
for the expected value of the degree as a function of fitness. The model alleviates two 
shortcomings that are present in the current models of growing multiplex networks: 
homogeneity of connections, and homogeneity of fitness. In this paper, we posit and 
analyze a growth model that is heterogeneous in both senses. 


keywords: multiplex networks, intrinsic fitness, grow¬ 
ing networks, preferential attachment, Stirling numbers. 


I. INTRODUCTION 

Multiplex networks are mathematical tools for mod¬ 
eling systems with multiple types of interaction. The 
system is conceptualized as being comprised of multiple 
layers, each hosting a distinct type of link (which corre¬ 
sponds to a type of interaction) between nodes. The set 
of nodes are the same for all layers. Many real systems 
have been modeled under the multiplex framework, such 
as citation networks oia, online social media [3] , airline 
networks [4] , scientific collaboration networks [1] , urban 
transportation networks [3, and online games [5]. The¬ 
oretically, multiplex networks demonstrate how incorpo¬ 
rating additional dimensions and types of interaction to 
simple one-layer systems can change their dynamics, add 
new properties and alter existing ones. Diverse processes 
have been studied theoretically on multiplex networks. 
Examples include epidemics [Hi, pathogen-awareness 
interplay [5], percolation processes [ziiiniiii], evolution 
of cooperation [I1II31, diffusion processes m and social 
contagion m- For thorough reviews, see min]. 

In the present paper we focus on the problem of grow¬ 
ing multiplex networks with fitness. Previous studies 
on growing multiplex networks exhibit two main short¬ 
comings: homogeneity of growth, and homogeneity (or 
absence) of nodal fitness. In [18], a growing two-layer 
network is studied, and various attachment kernels are 
envisaged. The number of links established by each new¬ 
comer is considered to be the same for both layers (In 
the Supplemental Material of [T8| , the possibility of het¬ 
erogeneous growth rates is entertained in the asymptotic 
mean-field analysis of degrees of individual nodes within 
each layer. However, the effect of growth heterogeneity 


on the single-layer and inter-layer degree distributions 
remains unknown). Similarly, the model posited and 
thoroughly analyzed in [19] exhibits homogeneity in the 
sense that, for a given node, the expected degree is the 
same across layers. In real systems, the nature of the 
connections in different layers differ, since they pertain 
to distinct types of interaction. For example, in [6|, the 
interactions between the players of a massive online game 
is mapped onto six distinct layers, and their average de¬ 
grees are different. It would be plausible to devise a 
growth model which incorporates heterogeneity explic¬ 
itly. In [201, the problem of homogeneity is alleviated 
by considering heterogeneous link growth rates. In the 
present paper, we consider heterogeneous link growth 
rates. 

Another unrealistic assumption that is made by the 
previous studies on growing multiplex networks is that 
the probability for each existing node to receive links 
from incoming nodes only depends on their degrees. For 
example, if we consider the network of citations between 
scientific papers, this assumption would mean that the 
inherent quality and novelty of the papers have no role in 
the future number of citations that they would receive. 
That is, only fame drives scientific success, not quality: 
when scholars cite a paper, they only take into account 
the number of citations that an existing paper has. This 
is obviously not the case. Similarly, consider the case of 
online social networks such as Twitter, Instagram, Pin- 
terest, Google-f and Tumblr. In all these networks, each 
user can ‘follow’ other users. Assuming that links are 
established only based on existing degrees—and not in¬ 
corporating any intrinsic fitness for the nodes—would be 
synonymous with disregarding the role of the quality of 
the content produced by each user on her/his popular¬ 
ity. True that after a user becomes famous, the fame 
on its own contributes to further accumulation of fol¬ 
lowers (which is the rationale behind all preferential at- 
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tachment models), but quality also has an undeniable 
role—especially, at the initial stages of the lifetime of 
each node (user). This motivates us to consider intrinsic 
fitness for nodes. In pTl425) . intrinsic fitness is envisaged 
in the case of single-layer networks. To our knowledge, 
no fitness-based model on multiplex networks exists in 
the literature. 


We consider a growing directed multiplex network that 
comprises two layers. Each node is assigned an intrinsic 
fitness, which models its quality. The fitness of a node 
never changes. Each node belongs to two layers: a merit 
layer and a fame layer. In the former, fitness values 
are the sole drivers of the growth mechanism. In the 
fame layer, attachment is preferential, that is, the prob¬ 
ability that a node receives a link from a newcomer is 
proportional to the total degree of that node, i.e., the 
sum of its degrees in both layers. For example, in the 
case of citation networks, the interpretation of the model 
is as follows. Two distinct types of citations can be dis¬ 
cerned. The first type—the meritocratic type—is when a 
scholar reads a paper, and cites it because of its content 
(a citation which would be given regardless of the num¬ 
ber of citations that paper already has). Another type 
of citation is what we call fame-driven. A paper can 
become trendy, or well-known in some literature (par¬ 
ticularly true for seminal papers which initiate a new 
subfield), and many citations that it receives would be 
solely due to its fame—i.e., current number of citations, 
which itself is the total of meritocratic and fame-based 
citations. For example, after a seminal paper initiates or 
revives a scientific domain, after the domain passes its 
inchoate stages, many papers will be remote from those 
seminal papers, but will still cite it because those pa¬ 
pers are famous, not because their content is being used 
(even tangentially) in the new paper being published. It 
is imperative to note that to a scholar who wants to cite 
an existing paper, quality is latent. That is, only the 
total number of links is observed; fame-based and mer¬ 
itocratic citations are not distinguishable for the new 
incoming node (the new paper). What is observable is 
the collapsed network, in which the links are aggregated 
into one layer. 


Our model emulates the said merit-fame interplay. We 
focus on the interlayer joint distribution of degrees and 
fitness. We find P{k,£,d), which is the (asymptotic) 
fraction of nodes with fitness 9 who have degree k in 
the merit layer and degree i in the fame layer. This is 
presented in Equation ( [T8| ). We also find P{q,6), which 
is the fraction of nodes with fitness 9 whose total degree 
is q. This is given in Equation (22). The results depend 
on the distribution of fitness values, as well as the initial 
number of links that each new node emanates in each of 
the layers. We also find the conditional expected total 
degree of nodes. That is, for a given fitness value, we 
find the expected number of total links. This result is 
presented in Equation (34). 

The rest of the paper is organized as follows. After 


introducing notation and terminology, we describe the 
growth mechanism quantitatively. We then undertake 
the rate equation approach to quantify the evolution of 
P{k, £, 9) as a function of time. We then focus on the 
steady-state, when transients vanish, and solve the re¬ 
sulting equations. We then obtain P{q, 9) through a 
straightforward transformation, and then use it to find 
the conditional expected value of total degree. 


II. NOTATION AND TERMINOLOGY 

The network is directed, and we use the terms degree 
and in-degree interchangeably. Form node x, the fitness 
value is denoted by 9^. The probability distribution of 
fitness values is denoted by p{9). The layer-I degree of 
node X is denoted by and the layer-2 degree of node 
X is denoted by ix- The total number of links of node 
X is denoted by qx, that is, qx = (-x + kx- If a quantity 
depends on time, we will explicitly mention it. If time 
dependence is not mentioned, the steady-state value of 
the quantity is meant. For example, kxft) is the degree 
of node x at time t, and kx is the degree of node x in the 
steady state, that is, in the limit as t —> oo. 


III. MODEL 

The system initially comprises N(fS) nodes, each with 
two types of links. The links are assumed to be estab¬ 
lished on two separate layers, layer I (the merit layer) 
and layer 2 (the fame layer). Let us emphasize that lay¬ 
ers embody the set of nodes, but the sets of links differ. 
Suppose that there are Li (0) links in the first layer and 
^ 2 ( 0 ) links in the second layer at the outset. 

The network grows by the successive addition of new 
nodes. Time increments in discreet steps, and at each 
timestep one new node is added to the network. Each 
incoming node establishes /3i layer-I links and /32 layer-2 
links to the existing nodes. 

Upon being born, the fitness of an incoming node is 
drawn from p(9) and stays the same thereafter. The 
mean value of the fitness distribution, that is, the ex¬ 
pected value of the fitness of incoming nodes, is denoted 
by p.. 

In the first layer, the probability of receiving a link 
from an incoming node for node x is proportional to 
9x: where 9x is the fitness of node x. In the second 
layer, the probability of node x receiving a link from the 
newcomer is proportional to kx £x- The probability of 
receiving a link in layer I can be written as ^ . The 

sum in the denominator can be computed at time t as 
follows. If the sum of the fitness values of the nodes at 
time t = 0 is 0, then as time progresses, the sum of the 
fitness values of nodes converges to Q-\-pt, where p is the 
mean of the fitness distribution. Since we will eventually 
limit the analysis to the steady state, the error of this 
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FIG. 1: Schematic illustration of the growth mechanism. A set 
of 7 nodes exist in the network, and node 8 is being appended. In 
the merit layer, only the fitness values drive nodes’ chances of 
receiving links from node 8 (who establishes /3i links in this 
layer). In the fame layer, the probability of receiving a link is 
proportional to total degree. For example, node 1 will have the 
highest chance of receiving a link, because it has the greatest 
total degree (which is equal to six). Conversely, node 7 has the 
smallest total degree (which is zero), and will not get a link from 
node 8 in the fame layer. 


approximation vanishes. For layer 2, the probability of 
receiving a link for node x is equal to y The 

sum in the denominator at any time equals the total 
number of links in both layers. 


IV. JOINT INTERLAYER DISTRIBUTION OF 
DEGREES AND FITNESS 


At time t, upon the addition of the new node, cer¬ 
tain events can change the value of Nt{k, £, 9). If a node 
of (fc — 1,£,9) receives a layer-1 link, then it becomes 
a {k,£,9) node, and Nt{k,i,9) increments consequently. 
Similarly, if a node oi {k,£ — 1,9) receives a layer-2 link, 
then it becomes a {k,£,9) node, and Nt{k,£,9) incre¬ 
ments consequently. On the other hand, if a node is 
already a (fc, £, 9) node and it receives a link in either 
layer, it will no longer be a {k,£,9) node, and Nt{k,£,9) 
decrements consequently. Also, note that the layer-1 de¬ 
gree and layer-2 degree of each incoming node is 0 upon 
introduction, and such a node has fitness 9 with prob¬ 


ability p{9). The following rate equation summarizes 
these events with their respective probabilities of occur¬ 
rence, where E{-} denotes expected value: 

E{Nt+i{k,9,£)} = Nt{k,9,£) 

, ^ 9Nt{k-l,£,9)-9Nt{k,9,£) 

+ - oTJd - 

{k + £- l)Ntik, £-l,9)-ik + £)Nt{k, 9, £) 
Li(0) + L2(0) + (/3i+/?2)t 
+ SkoSiop{9). ( 1 ) 

Hereinafter, we drop the expected value operator, 
and all the P{k,£,9) numbers denote expected values. 
Using the relation Nt{k,9,£) = (A^(0) + t)Pt{k,9,£), we 
can rewrite Q to quantify the evolution of Pt{k,9,£) as 
follows: 


[iV(0) + t] [Pt+i{k, 9, £) - Pt{k, 9, £)] + Pt+i(A:, 0, £) = 

^ 9Nt{k-l,£,9)-9Nt{k,9,£) 

Pi -o , - 


+ ^2 


{k + £-l)Nt{k,£-l,9) - {k + £)Nt{k,9,£) 

Ti(0) + L2(0) + (/?i+/32)t 


+ SkoSmp{S)- 


(2) 


Note that negative k or £ does not have a physical 
meaning, so k and £ in equation ([^, as well as every 
equation henceforth, are nonnegative integers. 

Now we focus on the steady state, where by definition, 
the values of Pt{k, 9, £) reach horizontal asymptotes and 
their variations vanish. Also we note that in the limit as 
t —>■ oo, we have 

lini ^ 

t-s-oo (d + pt a 

i: o _ ^(0) + ^ ^ 

t->-oo Li(0) -b ^ 2 ( 0 ) -f (/3l -f/32)t Pi + P 2 

(3) 

Using these limits, we can rewrite ([^ for the steady 
state as follows 

P{k-l,£,9)- P{k,£,9) 

? - l)P(fc, £-l,9)-{k + £)P{k, £, 9) 

+ SkoSeop{9). (4) 


P{k,£,9) = 

9- 


+ 


P 2 


Pi + P2 


{k + 


This can be rearranged and recast as 


P(fc,£,0) 



C-P) 


+ k + £ 


Pl{Pl+P2)0 p ^ 

P 2 P 

+ ^^^^Sko6iopi9). 


1,9) 

(5) 
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Dividing both sides by the factor on the left hand side, 
this transforms into 


P{kJ,B) = 

13iWi+P 2)0 \ Pik-1,£,0) 


+ 


+ 


{k + £-l)P{k,£-l,e) 


1 £i£ , , _ 

V + M M 


(i + ^) (^)+o + o' 


( 6 ) 


Hereinafter, for brevity of notation, we denote 
by A, and we denote ^1 + by Gg. Thus the differ¬ 
ence equation we need to solve takes the following form 

AOPi P{k-l,£,0) 


P{k,£,0)=- 


AGg + k + £ 

{k + £ - l)P{k,£ - 1,0) 6koSeop{0) 


AGg + k + £ 


Gg 


(7) 


. Let us define 


,p(k.e, 9) fe' (, 9). (8) 




It is easy to verify the following relations using the 
properties of the Gamma function: 

'0(fc — 1, £, 0) 


A0l3i P{k-1,£,0) 
p AGg + k + £ (^_^yY(AGg + k + £ + l) 

P{k, £ — 1,0) i’ik, £ — 1,0) 

AGg + k + £ 


(^) r(HG, + fc + £+l) 


(9) 


We substitute the first two terms on the right 
hand side of Q with the expressions given 
in Then we multiply both sides by the factor 

(aI^) + k + £+1). We arrive at the following 

difference equation: 

£, 0) =ij{k -l,£,0) + ik + £- l)^lj{k, £-1,0) 

+ ASkoS£op{0)T{AGg). ( 10 ) 

Without loss of generality, we can take k,k + £ to he 
the arguments of the function instead of k, £. Let us 
define the new auxiliary function: 

<l^l+\0) = Hk,£,0)- 


(Note that fc -I- £ is an upper index, not a power.) We 
can readily rewrite (10) in terms of (j). The difference 
equation reads 

£^l+\0) =0^+^l(0) + ik + £-l)4+^-\0) 

+ ASkoS£op{0)r{AGg). ( 12 ) 

Note that the last term on the right hand side is merely 
the boundary condition at {k,£} = {0,0}, and vanishes 
for any other combination oi k,£. As our last change of 
variables, let us denote k + £ — 1 by n. Dropping the 0 
argument for notational brevity, we can rewrite (12) as 


a function of n, k (without the term which dictates the 
boundary condition at {k,£} = {0,0}) in the following 
form: 




= +nq 


(13) 


This is the recurrence relation which defines the un¬ 
signed Stirling numbers of the first kind. We denote the 
Stirling numbers by [^] in this paper. Incorporating the 
initial conditions, the solution to (131 is 


X AT{AGg)p{0). 


Replacing n with k + £ — 1, we have: 


(14) 




'k + £-l 
k 


X AT{AGg)p{0). (15) 


Comparing this with we find the solution for 'll: 
to be as follows: 


^{k,£,0) = 


'k + £ 
k 


AT{AGe)p{0). 


(16) 


This readily yields P{k,£,0), using Q. We get 


P{k,£,0) = 


k + £ 
k 


Ar{AGg)pi0) 


r(AGg -\- k -\- £ -\-1) 

(17) 


Plugging in the explicit expressions for A and Gg, we 
arrive at the final solution: 


P{k,£,0) = 


'k + £ 
k 


Pl+ P2\ f (/?! + /?2)4 0 
1^2 ) \ h P 


r 

01+02 ( 1 1 /SiS} 

[02 V 


r 

Ai+02 

02 

{}+^]!)+k + £+l 


p{0). 


( 11 ) 


(18) 



















































V. COLLAPSED JOINT DISTRIBUTION OF 
DEGREE AND FITNESS 
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In real settings, the total number of links received by 
nodes are observed. For example, in the network of ci¬ 
tations between scientific papers, what is observed and 
documented is fc that is, the total number of links 
(citations) received by papers. One cannot observe the 
number of citations that a paper receives purely based 
on its merit (meritocratic attachment), or the number 
of citations it receives due to its popularity (fame-driven 
attachment). This motivates us to derive the distribu¬ 
tion of the total number of links, that is, k + i. Let us 
denote it by q. The joint distribution of q, k is simply 
P{k, q — k, 6). If we sum over all possible values of k, we 
get: 


P{q,e) = 


(ft+ft)r 


P 1 +P 2 

/32 


('+¥) 


'[* 5 ? (l + ¥)+« + l 

q 

E 


Pie) 




iPi + ^ 2)^1 e 


132 


(19) 


From the properties of the unsigned Stirling numbers 
of the first kind, the sum on the right hand side of (191 
can be readily evaluated. We have: 


P{q:0) = 


(ft^)r 


01+132 


1 + 


^) 


(/3i+fe)/3i ^ 


(92 


Pie) 


{Pi+MPi ^ 


02 




( 20 ) 


This is depicted in Figure 

If we use the generalization of binomial coefficients to 
non-integers, we can express ( |20| more concisely. Let us 
use the following notation for binomial coefficients: 


F(a-bl) 


F(6+l)r(a-6 + l)’ 


( 21 ) 


where a, b need not be integers. Using this notation, we 
can express (20) equivalently as follows: 


Piq,0) = 


Pi + P 2 

Pi + 2/32 


Pi+I32 

P 2 


(i + “)-i 


I3i+P2 


01+132 (1 , 

02 fJ. ) 

01+02 I 1 
02 


-Pie). ( 22 ) 


Now let us investigate the asymptotic behavior of 
Piq,e) for large values of q. From (20), we observe that 



q 0 


FIG. 2: Depiction of logarithm of inverse of P{q, 6). We have 
considered an exponentially decaying fitness distribution for 
illustrative purposes: p{9) = 0.9^/10. We have also set /3i = 2 
and ^2 = 4. The logarithm is taken for smoothing purposes, due 
to the rapid plummet of the function P{q, 9) in the q domain. 



FIG. 3: The log-log plot of P{q, 9) for three values of 9 for 
exponential fitness distribut ion p{9) = 0.9^/10, with pi = S and 
p 2 = 4. As we expect from l |24[ |, the curves are parallel for large 
values of q. It is visibly clear that the curves are also parallel to 
the asymptotic curve as indicated by \2A\ . 


only two terms have q. Using the Stirling approximation, 
we have 


PiqiO) 


iq + AGg A-\- l)<l+^Ge+A-\-3^-{q+AGe+A+l) ' 

(23) 


So we arrive at the following asymptotic relation 


P((?,6») - ye. (24) 

Note that the exponent does not depend on 0. 
This means that the degree distribution of the sub¬ 
populations with any fitness value follows the same ex¬ 
ponent. In other words, the rate at which the degree dis¬ 
tribution vanishes is the same for all fitness values. The 
relative chances of different nodes attaining extremely 
large degrees depend only on their fitness values, and 
not the degree itself, because if we divide the respective 
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probabilities, only the fitness-dependent multiplicative 
factors would determine the ratio, as the q-dependent 
parts cancel out. This is illustrated in Figure An¬ 
other implication of (24) is that the total degree distri¬ 
bution of the network, i.e. P{q), has a power-law tail 
with exponent 2 ^. 


VI. EXPECTED DEGREE DISTRIBUTION AS 
A FUNCTION OF FITNESS 

It is straightforward to compute the expected value of 
the degree distribution (p0|. We have 


{q)s = '^qPiq\0) = 


q=0 


q=0 


P{q,0) 

^ p{0) 


V J 


' Pl+h ( 

02 V 


1 + 


¥) 


■p ( (/^l+fe)/^! ^ 

‘ fi 


X ^<7- 

9=0 r 


Ylq+(Jl±^ 


q+l+^ 


1 + 




We now perform the following summation: 


„ def 

s = 


Ylq+iJl±^ 


0 


=0 r 


9=0 


1 I 01+02 
^ + 02 




(25) 


(26) 


In Appendix we prove the following identity for 
general real positive numbers y, x: 


E 

9=0 


r(g + y) ^_ 

r((3'-I-a;-I-y) (a: - l)r(y - 1-|-cc)' 


r( 2 /) 


If we use a: -I- I instead of a: in (27), we get 

r(g + y) 


E 

9=0 


_^ r(7/) 

r(y-I-a;-I-y-I-1) a;r(y-|-a:)' 


(27) 


(28) 


Now note that, using the basic properties of the 
Gamma function, we can rewrite (27) equivalently as 
follows 


E 

9=0 


{q + X + y)r{q + y) ^ _ 

r(y -I-a:-I- y -I-1) (a: - l)r(y - 1 -|- a;) 


r(y) 


■ (29) 


Expanding the left hand side, we have 

gr(g + y) ^ r(g + y) 


E 


^ r(g + a; + y-h 1) 

^ r(y) 

{x- l)r(y- 1-ha;) 


(a; -h y) E 


9^ ^(9 + ® + y + 1) 
(30) 


Combining this with (28), we arrive at 
gr(g -h y) 


E 

9=0 


ny) 


r(y-I-a:-I-y-I-1) (a; - l)r(y - 1-|-a;) 


- (a; -h y) 


r(y) 


:r(y- 


(31) 


Using the basic properties of the Gamma function and 
algebraic simplifications, we can express this in the fol¬ 
lowing form: 


V 9r(g + y) ^ r(y) y 

Yq r(g + a: -I- y + 1) F (x -|- y - 1) a: (a - 1) (x -|- y - 1)' 

(32) 


This has the same form as (26). 
tity ([^ with y = and x = 

S as follows: 


We can use iden- 
1° calculate 


-p { (01+02)01 e\ 

^ V 02 fl) 

p / ft I (01+02)01 
^ \02 ^ 02 ll) 

( ( 01 + 02)01 e\ 

V 02 + 

i ± _ 

( 01+02 \ f 01 \ ( 01 , (01+02)01 9 A ■ 
\ 02 ) {02 ) {02^ 02 +) 


Plugging this into (25), we get 


(33) 


{q)e = {Pi + ^ 2 )-- 


(34) 


This is a linear relationship (see Figure]^. If we take 
the average degree over all nodes, we need to sum up (341 


over all possible values of 9. In the numerator, y is 
created, which cancels out the y in the denominator and 
we get 


{q) = {Pi+P2)- (35) 

We know this result is true, because by construction, 
the total number of links created in the system (which 
is always equal to the sum of in-degrees of all nodes) is 
{Pi + P 2 )t at large times (when the effects of the initial 
conditions vanish) , and the total number of nodes is t, 
which means that their ratio (which yields the average 
degree) is equal to Pi + P 2 - 


VII. SUMMARY AND OPEN PROBLEMS 

This paper extends the literature of multiplex net¬ 
works by introducing a simple model which incorporates 
intrinsic fitness and preferential attachment. The merit 
layer is latent, yet drives the growth mechanism and the 
degree dynamics. We obtained closed-form expressions 
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FIG. 4: As | |34[ | predicts, the expected values of g is a linear 
function of 6, with slope and zero intercept. 


r(g + y) r(y) 

^ T{q + x + y) (x- l)r(y- 1 + a;)' 


The definition of the Beta function for positive real 
values x,y is 


B{a,b)= [ = 

Jo 


r(Q)r(b) 
r(a + b) ' 


(A2) 


We can rewrite the summand of (All as follows: 


for the joint interlayer distribution of degrees and fitness, 
as well as that of the total degrees. We observed that 
the expected value of the total degree linearly increases 
with fitness. 

An immediate generalization of the present problem 
would be its extension to an arbitrary number of lay¬ 
ers. Also, we have disregarded the temporal dynamics 
of the system and its transients in favor of the steady 
state. This loses valuable information about the tran¬ 
sient state and the effects of initial conditions on the 
evolution of the network. In other words, for a given ini¬ 
tial network (not necessarily small), one can study the 
evolution of the system in arbitrary time regimes, and 
investigate how the properties of the initial network af¬ 
fect the equilibration of the system, and how they affect 
the asymptotic properties of the network. 

Another immediate step to augment the present model 
is to devise statistical recipes for inference. Since fitness 
is a latent variable and only q can be observed, one can 
use ( [20| ) (or its time-dependent version) to devise max¬ 
imum likelihood techniques to infer the distribution of 
fitness (merit) of scientific publications, blog posts, etc., 
by observing the distribution (or evolution) of degrees. 


Appendix A: Proof of Identity (271 


r(g + i/) 


Jo 


r(g 


y) 


r(x) 


(A3) 


We now have 


E 


r(g + 2/) 


1 


OO „1 

/ dt{i - 

_n Jo 


^^^T{q + x + y) T{x)—j^ 

Jo ^0 

r(a:) Jo 1-t 

i [X) 1 (x) I (x -I- 2/ — 1) 

1 r(x-i)r(y) 


(x - i)r(x -1) r(x -I- ?/ - 1) 

ny) 

(x- l)r(y- l+xY 


(A4) 


We need to prove the following identity 


which concludes the proof. 
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